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Abstract. We investigate the accuracy of mass determinations Mbh of supermassive black holes in 
galaxies using dynamical models of the stellar kinematics. We compare 10 of our Mbh measure- 
ments, using integral-field OASIS kinematics, to published values. For a sample of 25 galaxies we 
confront our new Mbh derived using two modeling methods on the same OAS I S data. 
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INTRODUCTION 

Scaling relations between global parameters of galaxies and their central supermassive 
black holes (BHs) [1, 2, 3, 4, 5, 6] are central to our current understanding of how galax- 
ies assemble in the hierarchical merging paradigm [7, 8]. Thousands of publications 
made use of a small collection of ~ 50 BH mass determinations Mbh accumulated from 
a large combined effort over the past 15 years (see compilations in [9, 10, 11]). Given 
the relevance of Mbh measurements, and the still small numbers available, it is essential 
to progress with the observations and models of new galaxies to populate the scaling 
relations. Together with the relations themselves, their intrinsic scatter is also important. 
It is used to test models of the joint evolution of BHs and galaxies and is required to 
estimate the mass function of supermassive black holes in the Universe [12, 13, 14]. For 
the latter one needs a reliable estimate of the measurements errors. Here we present an 
update on our efforts to address these issues: (i) Populate the BH scaling relations and 
(ii) test the accuracy of the Mbh measurements. 



MODELS AND RESULTS 

Our models are constrained by high spatial resolution (but seeing-limited) integral- 
field observations of the stellar kinematics obtained with OASIS [15], in combina- 
tion with large field SAURON kinematics [16]. We adopt the Multi-Gaussian Expan- 
sion (MGE) parametrization of the surface brightness [17]. We apply two independent 
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FIGURE 1. Data-model comparison for the Schwarzschild model. The bi-symmetrized OASIS stellar 
kinematics [15] of NGC 4473 (top panels) is compared to the best-fitting axisymmetric Schwarzschild 
model (bottom panels). From left to right the plots show the mean velocity V, the velocity dispersion a 
and the Gauss-Hermite moments hj-h^- Contours of the galaxy surface brightness are overlaid. 



axisymmetric modeling methods to fit the same stellar kinematics and MGEs: (i) An 
orbit-superposition Schwarzschild's [18] approach as implemented in [19] and (ii) the 
Jeans Anisotropic MGE (JAM) method [20]. The availability of good OASIS kinemat- 
ics, Hubble Space Telescope (HST) imaging and the exclusion of objects clearly incon- 
sistent with the assumed axisymmetry lead to a sample of 25 early-type galaxies. 

An example of the results obtained with our two modeling approaches on NGC 4473 
is shown in Fig. 1, 2 and 3. This galaxy is an interesting case as it shows no peak 
in the central V rsm (or a) and one may wonder whether a BH is really present and 
what observable is actually constraining its mass. It was also modeled by [21] who 
derived M BH = (l^l^) x 10 8 M (la errors), using HST/STIS kinematics. The value 
is consistent with our OASIS determinations, but it has larger errors. Even though the 
radius of the BH sphere of influence is 7?bh ~ 0.24 /; while our seeing PSF has a 3.3 x 
larger FWHM rs 0.80", the BH can be reliably measured. Fig. 2 illustrates the fact that 
the BH detection in this object comes from the lack of the clear drop in V rsm that one 
would expect to see if there was no BH, with our spatial resolution. The integral-field 
data are essential to constrain the anisotropy and provide a clear BH detection. 

In Fig. 4 we confront our Mbh determinations, using Schwarzschild's method and 
OASIS data, to 10 measurements from the compilation of [11] using the same method 
but HST/STIS kinematics. There is a reassuring agreement between the two completely 
independent determinations. However the differences between the two sets of values 
would imply a la error on Mbh of ~ 0.20 dex, which is significantly larger than 
the typical formal errors. This is likely due to unaccounted systematics in the data 
and differences in the modeling implementations. We measure a similar scatter when 
comparing all our 25 Mbh values with Schwarzschild's method against the ones obtained 
with JAM. As the former method is more general than the latter, this scatter gives an 
upper limit to the errors in our OASIS determinations, however it does not take into 
account systematics in the data and in the modeling assumptions. Below Mbh ^ 10 8 
M most of our OASIS measurements are upper limits. The low Mbh regime must be 
explored with observations using adaptive optics as illustrated in this volume by [22, 23] . 




FIGURE 2. Data-model comparison for the JAM models. The first three rows show the Vrms of the JAM 
models for different M B h (in units of 10 8 M©) and anisotropy f} z = l — <5 2 /o\ (above each plot), with 
M/L optimized to best fit the OAS I S data. The bottom-left panel is the bi-symmetrized V rms = y/V 2 + a 2 
observed with OASIS. The best-fit can be done by eye and it agrees with the A.% 2 minimum of Fig. 3. 
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FIGURE 3. Left Panel: Contours of A^ 2 describing the agreement between the OAS I S + SAURON data 
and the Schwarzschild dynamical models, as a function of the Mbh and M/L. Right Panel: Adaptively- 
sampled contours of A^ 2 between the OASIS data and the JAM models, as a function of Mbh and the 
anisotropy j5 z . For each model the M/L was optimized to best fit the OASIS data. The three lowest 
contours correspond to the 1, 2, and 3(7 (thick line) confidence levels, marginalized for one parameter. 
The black dots show the models that were actually computed, while the contours were interpolated with a 
thin-plate spline. 
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FIGURE 4. Left Panel: Comparison between the Mbh measured with our Schwarzschild models of the 
OASIS integral-field kinematics and published Mbh from [11] using HST/STIS kinematics. Upper limits 
in our models are indicated by an arrow. Right Panel: Comparison between our Mbh from Schwarzschild 
and from JAM models of the same OAS I S data. Upper limits are indicated by star symbols. 
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